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ABSTRACT 

We explore the behaviour of [C n] A157.74 fim forbidden fine-structure line observed in 
a sample of 28 galaxies selected from ~ 50 deg 2 of the H- ATLAS survey. The sample 
is restricted to galaxies with flux densities higher than S i 60 ^m > 150 mJy and opti¬ 
cal spectra from the GAMA survey at 0.02 < z < 0.2. Far-IR spectra centred on this 
redshifted line were taken with the PACS instrument on-board the Herschel Space Ob¬ 
servatory. The galaxies span 10 < log(LiR/L Q ) < 12 (where Lir = Lir[ 8 — 1000/im]) 
and 7.3 < log(L[cn]/L 0 ) < 9.3, covering a variety of optical galaxy morphologies. 
The sample exhibits the so-called [Cn] deficit at high IR luminosities, i.e. F[cii]/£ir 
( hereafter [C ii] /IR) decreases at high Lir. We find significant differences between 
those galaxies presenting [Cn]/IR> 2.5 x 10 -3 with respect to those showing lower 
ratios. In particular, those with high ratios tend to have: (1) Lir < 10 11 L 0 ; (2) cold 
dust temperatures, T < 30 K; (3) disk-like morphologies in r-band images; (4) a WISE 
colour 0.5 < Si 2 ^m/S 22 gm ^ 1.0; (5) low surface brightness Sir « 10 8-9 L e kpc -2 , 
(6) and specific star-formation rates of sSFR^ 0.05 — 3Gyr _1 . We suggest that the 
strength of the far-UV radiation fields (<Go>) is main parameter responsible for 
controlling the [Cii]/IR ratio. It is possible that relatively high <Go> creates a pos¬ 
itively charged dust grain distribution, impeding an efficient photo-electric extraction 
of electrons from these grains to then collisionally excite carbon atoms. Within the 
brighter IR population, 11 < log(Lni/L 0 ) < 12, the low [Cii]/IR ratio is unlikely 
to be modified by [C n] self absorption or controlled by the presence of a moderately 
© 2015 RAS luminous AGN (identified via the BPT diagram). 

Key words: Galaxies: starburst; ISM lines and bands; ISM evolution; Infrared: ISM 
resolved and unresolved sources as a function of wavelength; 


1 INTRODUCTION 

Understanding the chemical and physical evolution of a 
galaxy is far from trivial. Newly born stars consume 
and process the available gas, whilst heating the inter¬ 
stellar medium (ISM), and supernovae enrich the envi¬ 
ronment with heavy elements, contributing to potentially 
complex feedback processes. A good description of ISM 
physics under the influence of stellar radiation fields was 
achi eved using photo-dissociatio n region (PDR) modelling 
(e.g. iTielens fc Hollenbachll 19851 ). These models can explain 
the origin of most of the dense ISM emission from star¬ 
forming galaxies, including the major cooling fine-structure 
lines of carbon/nitrogen/oxygen (C/N/O) and the underly¬ 
ing Infrared (IR) continuum emission produced by interstel¬ 
lar dust. 

At IR wavelengths, the most prominent emission line is 
[C n] A157.74 /im CP 3/2 ~^ 2 Pi/ 2 i Pui/k = 92 K), which car¬ 
ries ~ 0.1 — 1 per cen t of the bolometric pow er emitted by 
star-forming galaxies (IStacev et alJll99ll . l2010h , with a [Cll] 
flux typically lOOOx that of CO (J = 1-0) at 115 GHz. The 
low ionisation potential, 11.26 eV, makes [C11] a key partic¬ 
ipant in cooling the warm and diffuse media, converting it 
into c old and dense clouds that can then collapse to form 
stars dDalgarno &; McCravl 119721 ). As a fine-structure line, 
[C 11] is an excellent tracer of all the different stages of evolu¬ 
tion of the ISM: it can be excited by collisions with electrons 
in the warm ionised medium; HI in the warm or cold diffuse 
media; and H2 in the warm and dense molecular gas. Its 
intensity is sensitive to the column density, th e volume den ¬ 
sity a nd the kinetic temperature of the ISM (IPineda et al.l 

120R . _ 

In the plane of the Milky Wav. lPineda et ahl 1120131 1 show 
that the [C n] line emission emerges predominantly at Galac- 
tocentric distances bet ween 4 and 1 0 k pc. Co nsidering a scale 
height for the Galaxy, IPineda et al.l (12014 ) finds that the 
ISM components that contribute to the [C n] luminosity of 
the Galaxy have roughly comparable contributions: 30 per 
cent comes from dense PDRs, ~ 20 per cent from ionised 
gas, ~ 25 per cent from diffuse atomic gas and ~ 25 per cent 
from CO-dark H2 (at the surface of molecular clouds where 
carbo n is not in the form o f CO). In the local spiral galaxy 
M 33, 1 Kramer et al.l (|2013l ) suggest that the [C11] emission 
related to neutral gas corresponds to ~ 15 % of the total 
within 2 kpc from the galaxy centre, while this percentage 
seems to increment up to ~ 40 d= 20 % in the outer part of 
the disk (between 2 and 7 kpc). 

In distant galaxies, strong observational limitations im¬ 
pede a detailed characterisation of the different phases of 
the ISM. Typically detected in a single telescope beam, the 
ISM phases are all mixed together, hence [Cll] line detec¬ 
tions relate to averaged quantiti es of an ensemble of indi¬ 
vidual PDRs, ionised regions, etc. lMadden et all |l993) sug¬ 
gest that ~ 75 % of the [C 11] emission from the spiral galaxy 
NGC 6946 originates from cold neutra l hydrogen clouds, an d 
<40% from the diffuse galaxy disk ([Contursi et al.|[200l j. 
On the other hand, averaged ISM properties of luminous- 
and ultra-luminous- IR galaxies (LIRGs and ULIRGs) are 
found to vary considerably with respect to those observed in 
the Milky Way. LIRGs/ULIRGs present much higher star- 
formation rates and evidence mu ch larger amounts of io nised 
gas - up to 50 % of the total (e.g. lMalhotra et ai1l200ll ). The 


rest of their [C11] lumi nosi t y is expected to come from dense 
PDR-r elated ISM (e.g . Negishi et al.ll200il ). At the other ex¬ 
treme, I Cormier et al.l (l2012l j estimate that the low metallic- 
ity dwarf star-forming galaxy Haro 11 produces only 10% 
of its [C 11] emission in PDRs, probably because radiation 
fields penetrate deep into the ISM components. 


The examples shown above clearly evidence the in¬ 
tricate decomposition of the [C11] emission into the dif¬ 
ferent ISM phases. Great effort has been dedicated to 
find [C 11] correlations with various galaxy properties. Pre¬ 
vious studies have shown that the integrated [C11] line 
strength depends on the Poly cyclic Aromatic Hydrocar- 
bon (P AH) emission lines (e .g. iBakes &; Tielensl 1 1994t h IR 
colour (Mal hotra et al. |200llh the degree of active galac¬ 
tic nuclei (AGN; e.g. Stacey et al.l l2010l ) activity, among 
other various correlations. Those with powerful AGN show 
fainter [C ll] than p u re star-forming galax ies, at fixed Lir 
(iNegishi et al.l l200ll : ISargsvan et al.l 120121 ) . [C11] has also 
been used for dia gnostic purposes in galaxies at cosmologi¬ 
cal distances (e.g.Tlvison et alJl20ld : IValtchanov et al1l201ll : 
I George et al.l I2OI3I) and has somet i mes betrayed a gal axy ’ s 
redshift (e.g. ISwinbank et ~aJ~]l2Q12l : iGeorge et al.ll2013l ). In¬ 
deed, [C11] is a potentially unrivalled tracer of the total gas 
mass or the star- format ion rate or the dynamics in the most 
distant galaxies (Maiolino et alJ i 2005: D e Looze et alJl2Q14l : 


iDe Breuck et al.l 20141: 1 Herrera-Camus et al.l 20151). 


Early observations of local st ar-forming g a laxies with 
the Kuiper Airborne Observatory (Sta cey et al.l 1991) and 
the Infrared Space Observatory (ISO: [Malhotra et al.l 1997 ) 
showed that [Cii]/IR is roughly constant at 1-3 xlO -3 
(with a factor of three of scatter) for Lir < 10 11 L 0 

galaxies (it is higher in low-m etallicity environments - 
lMaddeiJl2000l : [Rubin et ahl 1 20091 ) . although it drops rapidly 
at higher Lir. This behaviour is usually referred to as 
the ‘[Cll] deficit’. [Cn]/IR also correlates strongly with 
60 £un /100 pm col our (a proxy for dust temperature, e.g. 
iDfaz-Santos et ahl [20131) and the IR/L-band luminosity ra¬ 
tio, Lir/Lr ( Malhotra et al.ll200ll ). Various explanations for 
the [C11] deficit have been offered: (1) in high radiation 
fields, photoionisation of dust grains might saturate, hence 
the energy of their ejected p hoto-electrons is red uced, re¬ 
ducing also the gas heating (iLuhman et al.l [20031 ) : (2) the 
deficit is due to an increase in the collisional de-excitation 
of the rC i i] t ransit ion due to an increase in gas density 
([Negishi et al.l 200 ll ) : (3) a significant fraction of the IR 
emission may arise from dust absorption of photons from old 
stella r populations which are not related directly to PDRs 
(e.g. iRowan-Robinson et al.l 12010 ): (4) a significant por¬ 
tion of the IR emission could emanate from ‘dust-bounded’ 
structures wit hin photoionised gas regions, but the [C ll] 
line does not (ILuhman et al.l 1 20031: lAbel et al.l 120091) : (5) 
the [C 11] lin e may be self-absorb ed or optically thick (e.g. 
NGC 6334 - Boi^kp_^ 3etdll995l ) in highly embed ded re¬ 
gions ([Fischer et al.ll2010l ) (see also lGerin et ahl l2015l b or (6) 
the level of IR emission in the most IR-luminous objects may 
be bo osted by AGN activity (ICurra n 2009; Sargsvan et al.l 

l2012h . 


Recent evidence, based on observations of local 
star-forming galaxies with the Herschel Space Observa- 

































































































tor $_ ] (|Pilbratt et al.l l2010l b points to different modes 
of star formation dri ven by the star- formation effi¬ 
ciency, Lir/Mh 2 (e.g. lYoung et al.l Il986l b regardless of 
the origin of the ionised or neu t ral phase of the ISM 
(|Gracia-Carpio et alJ 12011). iGracia-Carpio et al.l find that 
Lir/Mh 2 ^ 100 L© Mq 1 marks a point at which the average 
properties of the neutral and ionised gas are significantly 
different. Using a sample of powerful star-forming galaxies 
they propose a scenario in which highly compressed, more 
efficient star formation, creates largely enhanced ionisation 
parameters that manifest themselves in lower line to con¬ 
tinuum ratios. This value o f Lir/ Mh 9 is closely related to 
that at which iGenzel et al.l (|2Q10h and iDaddi et al.l (l2010h 
claim a transition to a more efficient star-formation mode, 
above the so-called ‘main sequence’ for star-forming galaxies 
([Elbaz et al ]IB)- 


Herschel was able to observe [Cli] line emission from 
many local galaxies - a legacy that will last into the 
era of SPICA (the Spac e Infrared Telescope f or Cosmol¬ 
ogy and Astrophysics - Nak agawa et al.l l2012h . Compre¬ 
hensive analyses of these observations are mandatory if 
we are to interpret ground-based observations of [C 11 ] to¬ 
wards high-redshift galaxies, e.g. with the Atacama Large 
Millimeter Array (ALMA), where the [Cill line is shifted 
into accessible atmospheric windows (e.g. | Maiolino et al.l 
200.1120091: IWalter et, al.l i2 00i iDe Breuck et, al.ll201ll. 12014 


Riechers et al J 120131 . 2014h . Indeed, recent studies of lu¬ 


minous, high-redshift galaxies have revealed a [Cn] be¬ 
haviour which contrasts with that seen locally. Instead of 
looking like powerful ULIRGs - at the same luminosity 
- they exhibit striking similarities in terms of [C il l /IR 
ratio s to normal local star-forming g alaxies (|lvison et al.l 
20ld: lH ailev-D unsheath et al.l 1201 Ol : IStacev et al.l l20ld : 
Valtchanov et al.l 201 il l. 


In this paper, we present Hersc hel Photodetector Ar - 
ray Camera and Spectrometer fPACS: IPoglitsch et al.ll2010h 
Integral Field Unit (IFU) detections of a sample of dusty 
galaxies. We populate the L[ C n\ versus Lir diagram with 
galaxies that lie at distances between 90Mpc (z ~ 0.02) 
and 1000 Mpc (z ~ 0.2). This study therefore bridges the 
gap between local galaxies and the growing body of data 
acquired for galaxies at high redshift. Our study performs 
a uniquely wide and detailed parameter space exploration 
that perfectly complements recent Herschel studies at low 
and inter mediate redsh i fts, e .g. the Herschel ULIRG Survey 
fHERUS: iFarrah et alJl2013h an d the Great Observator ies 
All-sky LIRG Survey (GOALS - iDlaz-Santos et al ]|20l3),a 
sample of LIRGs at z < 0.09. 

In what follows, we adopt a K ronpa initial mass function 
(IMF - iKroupa & Weidnerll2003h and a ACPM cosmology 
with Ho — 70kms -1 Mpc -1 , Om — 0.27 and — 0.73. 


2 SAMPLE SELECTION FOR PACS 
SPECTROSCOPY 

We make use of the internal pha se-1 data release (v2; 
iRigbv et al.ll201ll : ISmith et al.ll201ll . Valiante et al. in prep., 


1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 



Figure 1. IR luminosity versus redshift for the 28 galaxies that 
were observed with PACS for [Cli] spectroscopy (large red cir¬ 
cles) and the full parent sample of 84 sources (black empty 
circles). Our targets cover approximately two orders of mag¬ 
nitude in Lir. We overp lot the Revised Imperial IRAS -FSC 
Redshift Catalogue of IWang et al.l ( 2014 1 [blue crosses; consid¬ 
ering L (8 — 1000/i,m)/L© = 47rD 2 x Sir, where Sir = 1.8 x 
1 0 -14 (13.48 Si 2 nm + 5.16 S 25a .rn + 2.58 S 6 o M m + Sioo^m) w m- 2 ; 
ISanders Sz Mira bel 1996], and the observed H -ATLAS sources in 
the GAMA-09h field (green crosses). The dotted line shows the 
detection threshold used to select our targets, above which sources 
were randomly selected (any objects lying below the threshold 
are due to an improved SED fitting approach that was imple¬ 
mented post-selection, resulting in slightly lower Lir estimates). 
The detection threshold was applied to ensure a > 5cr line detec¬ 
tion within 10 min using the simplest (single scan, single pointing) 
PACS spectroscopic mode, i.e. S^cn] > 15.5 X 10 -18 W m -2 as¬ 
suming [Cii]/IR = 0.001. 


and Bourne et al. in prep) taken from the three equatorial 
fields of the ijersc/ze LAstroph ysic al Te rahertz Large Area 
Survey (i/-ATLAffi lEalcs ct alJ hold . We make use of 
both Spectral and Photometric Imaging Rec eiver (SPIRE; 
Griflfin et al.ll2010l) and PACS H -ATLAS data ([Pascale et al.l 
20 111 : llbar et al.l ‘ 2010h . From the total number of 109,231 
sources detected in these fields, we performed the follow¬ 
ing selection criteria for our targets: (1) a flux density 
threshold at Si 60 /xm > 150 mJy, where 160 fim is near the 
peak of the spectral energy distribution (SED) of a lo¬ 
cal star-forming galaxy; (2) targets without PACS 160-/xm 
(> 3cr) neighbours within 2arcmin from their centroids (to 
avoid problems when chopping; see §[3]); (3) an unambigu¬ 
ously i dent i fication in the Sloan Digita l Sky Survey (SDSS 
DR6 -lA delman -McCarthv et al.ll2008h (reliability > 0.8, 
ISmith et al.ll20lH " Bourne et al. in prep.); (4) galaxies need 
to be smaller than the PACS spectroscopic field of view, so 
restricted to sources with Petrosian SDSS radii smaller than 
15arcsec in the r-band (see Fig. [3]); (5) high-quality spec¬ 
troscopic redshifts from the Galaxy and Mass Assembly sur- 
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vey (OAMAFI - iDriver et, al.lEifjl. 1201 ill (z_QUAL^3). Note 
that GAMA combines spectra from SDSS with deeper spec¬ 
tra taken with the AAOmega fibre-fed spectrograph on the 
3.9-m Anglo-Australian Telescope; ( 6 ) a redshift between 
0.02 < z < 0.2 (median, 0.05), where the upper limit is im¬ 
posed by the point at which the [C n] emission is redshifted 
to the edge of the PACS spectrometer 160 /im band (which 
is known to leak), and the lower limit is imposed by the 
Petrosian r-band criterion. Applying all these criteria, we 
remain with a sample of 327 sources. 

For each held, we randomly selected galaxies to span 
a wide range of optical morphological types and IR lumi¬ 
nosities. In total, we selected a statistically significant sam¬ 
ple of 84 galaxies for a [C n] spectroscopic campaign us¬ 
ing Herschel- PACS during its second open time call. Galax¬ 
ies have photometric detections at 100, 160, 250, 350 and 
500 /im by Hersch el , including Wide -held Infrared Survey 
Explorer ( WISE - I Wright et al.ll2010h photometry, and ap¬ 
proximately half of them were also detected by IRAS at 
60 /im (using the Revised I mperial IRA S- FS C Redshift 
Catalogu<0, RiFSCZ_SHORT.v2: IWang et al1l2014 ). The com¬ 
bination of IRAS/ WISE/Herschel/SDSS photometry and 
GAMA spectroscopy allows accurate estimates for a wide 
range of physical parameters, including Lir (having a range 
of 10 10 < Lir/Z /0 < 10 12 ), dust temperature, stellar masses, 
emission line strengths, metallicities, etc. (see § 01 ). 


3 PACS IFU OBSERVATIONS AND DATA 
ANALYSIS 

Data were obtained between 2013 April 18-28, a few weeks 
before Herschel ran out of liquid helium, for the Herschel 
Open Time project, ot2_eibar_1 (P.I. E. Ibar). Of the par¬ 
ent sample of 84 galaxies, 28 were observed (Table [3]), all 
from only one of the equatorial if-ATLAS coverages - the 
GAMA-09 h held. The selection of these 28 was performed 
purely on the basis of scheduling efficiency, so this sample is 
representative of the parent sample, just smaller in number 
(see Fig. [1]). 

We targeted the redshifted [C n]-157.7 /im line emission 
using PACS spectroscopy in the hrst order/Rl hlter. PACS 
comprises an array of 5 x 5 spaxels, 9.4arcsec x 9.4arcsec 
each, covering ~ 2 /im of bandwidth covered. The instanta¬ 
neous held of view in a single pointing (pointed-mode) ob¬ 
servation is thus 47 arcsec x 47 arcsec, hence based on the 
optical size, the majority of the line emission falls on the 
central spaxel. 

The observations were made with a small chopping an¬ 
gle (1.5arcmin), as appropriate for our sources. Note that 
our selection criteria have ensured that no sources would lie 
in the chopped beam. The H- ATLAS survey targeted areas 
of low IR background, so that the background intensities 
never approached the detector saturation level. 

We retrieved the calibrated PACS level -2 data prod¬ 
ucts (processed with SPG vl2.1.0) using the Herschel User 
Interfac43 (HUI V6.0.4). We exported the hps3drb (blue) 


3 WWW.GAMA-SURVEY.ORG 

4 astro.ic.ac.uk/home/mrrobinson 

5 ARCHIVES.ESAC.ESA.INT/hSA/ui/hUI.JNLP 


Table 1 . Log of PACS spectroscopic observations. Each obser¬ 
vation lasted 394 s and was processed with SPG version 12.1.0. 
The table shows the IAU name, H- ATLAS’s nickname, GAMA 
identifier and the OBSID during the Herschel campaign. 


IAU name 

Nickname 

GAMA ID 

OBSID 

HATLAS J09:17:21.91+00:19:18.8 

G09.v2.117 

601323 

1342271048 

HATLAS J09:12:05.82+00:26:55.5 

G09.v2.42 

216401 

1342270763 

HATLAS J09:09:49.59+01:48:45.9 

G09.v2.26 

324842 

1342270761 

HATLAS J09:07:50.13+01:01:42.6 

G09.v2.107 

279387 

1342270762 

HATLAS J09:05:32.66+02:02:20.0 

G09.v2.58 

382362 

1342270757 

HATLAS J09:00:04.98+00:04:46.7 

G09.v2.55 

209807 

1342270755 

HATLAS J08:58:35.96+01:31:48.9 

G09.v2.76 

376679 

1342270657 

HATLAS J08:58:28.62+00:38:14.8 

G09.v2.80 

622694 

1342270756 

HATLAS J08:57:48.00+00:46:41.2 

G09.v2.48 

622662 

1342270656 

HATLAS J08:54:50.33+02:12:08.9 

G09.v2.38 

386720 

1342270658 

HATLAS J08:54:06.05+01:11:30.5 

G09.v2.137 

301346 

1342270655 

HATLAS J08:53:56.59+00:12:55.6 

G09.v2.170 

600026 

1342270649 

HATLAS J08:53:46.47+00:12:51.6 

G09.v2.45 

600024 

1342270648 

HATLAS J08:53:40.87+01:33:48.1 

G09.v2.103 

323855 

1342270654 

HATLAS J08:52:34.39+01:34:19.8 

G09.v2.87 

323772 

1342270653 

HATLAS J08:51:12.83+01:03:43.6 

G09.v2.235 

371789 

1342270651 

HATLAS J08:51:ll.48+01:30:06.9 

G09.v2.60 

376293 

1342270652 

HATLAS J08:49:07.15-00:51:40.2 

G09.v2.175 

3624571 

1342270647 

HATLAS J08:46:30.79+00:50:55.1 

G09.v2.90 

278475 

1342270650 

HATLAS J08:44:28.41+02:03:49.8 

G09.v2.52 

345754 

1342270371 

HATLAS J08:44:28.27+02:06:57.3 

G09.v2.77 

386263 

1342270370 

HATLAS J08:43:50.90+00:55:34.0 

G09.v2.102 

371334 

1342270372 

HATLAS J08:43:05.18+01:08:57.0 

G09.v2.167 

300757 

1342270373 

HATLAS J08:42:17.71+02:12:22.3 

G09.v2.232 

867786 

1342270365 

HATLAS J08:41:39.45+01:53:44.8 

G09.v2.299 

345647 

1342270364 

HATLAS J08:38:32.01+00:00:44.5 

G09.v2.lll 

208589 

1342270374 

HATLAS J08:37:45.33-00:51:42.3 

G09.v2.66 

3895257 

1342270362 

HATLAS J08:36:01.57+00:26:18.1 

G09.v2.23 

214184 

1342270363 


and hps3drr (red) data cubes, including the appended sig¬ 
nal, coverage, noise and wavelength index. These data cubes 
comprise 5x5 spaxels x n c h an rebinned spectral channels 
for each observed galaxy. At the observed frequencies, the 
effective spectral resolution is ~ 190-240 km s -1 , providing 
useful kinematical information for typical disk-like galaxies 
(see Table[3]). We use the Interactive Data Language (idl) to 
analyse these cubes along with the wavelength index array. 

First, we find the best flux density measurement by 
comparing the central spaxel line emission (plus an aperture 
correction ~ 0.4-0.5; a point-like estimate) with the summed 
over the whole IFU (extended estimate). We conclude that 
the best compromise in terms of signal-to-noise for the line 
flux density is to use the weighted sum (aided by the ap¬ 
pended instrumental noise cube) of the central 3x3 spaxels. 
The addition of the outer spaxels only introduced noise in 
the signal. We recover, on average, 10% more flux in the 
extended estimate than in the point-like one. 

[C n] line fluxes were measured via a Gaussian fit to 
the added spectra (see Fig. El, fitting simultaneously a lin¬ 
ear background slope and a Gaussian using the mpfitpeak 
routine within idl. To perform the fit, we removed channels 
suffering from higher noise in both ends of the spectrum. 

Using the spect roscopic redshifts, we calculate [ C n] lu¬ 
minosities following ISolomon &; Vanden Boutl (2005). In or¬ 
der to estimate uncertainties for the line parameters, we 
ran a Monte-Carlo realisation (lOOOx), randomly varying 
the signal per spectral channel using the instrumental er¬ 
ror cube. Based on these simulated data, we quote l-a un¬ 
certainties for the line measurements and 3 -a upper limits 
(based on intrinsic 400 km s -1 FWHM widths) in Table [3] 

We opted to receive data in order/B3 (blue filter at 
70/im) rather than order/G2 (green), as we already have 
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Figure 2. The flux densities versus observed wavelengths for all 28 galaxies presented in this work. From short to long wavelengths, 
SED data points correspond to broadband photometry from: WISE 12 and 22 /im, IRAS 60 /im, Herschel PACS 100, 160, and Herschel 
SPIRE 250, 350, and 500 fi m. Note that WISE 12 is not used to fit the SED. All photometric points have significance above 3-cr limit. 
The inset in each panel shows the [C il] spectra, flux density versus rest-frame wavelength, observed by the central spaxel of the PACS 
detector. Only two galaxies (G09.v2.87 and G09.v2.167) resulted with a [Cn] peak to noise line ratio lower than three. These two galaxies 
appear as upper limits in the following figures. 
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Figure 3. Postage-stamp images (P X P) for our targets taken from the SDSS r-band imaging. 



Figure 4. Left: [Cll] line luminosity versus bolometri c IR luminosity for the s amp le of sources pre sented in this wo rk. We compare ou r 
observations to previous (at similar redshifts), GOALS (iDfaz-Santos et ahll2013h and lSargsvan et a! 1 ll2012l)’s samples JPfa^Rjmtos et al.l ’s 
data have been scaled by a 1.75 factor to convert from far-IR to IR luminosities (see Appendix E from lHerrera-Camus et . Con- 

tinuum and dashed red lines show the best linear fit (excluding [Cll] undetections) logio(L[cn] /±o) = (0.721 ±0.003) X log 10 (LiR/L©)± 
(0.53 ± 0.04) and the 0.2 dex scatter, respectively. Diagonal lines show the 1.0, 0.1 and 0.01 per cent fraction of Lir. Right: [Cll]/LiRas 
a function of Lir, revealing the well-known [Cll] deficit. Colour coded points are the same in both figures. 


100 -/im continuum photometry from the if-ATLAS imaging 
survey. This PACS blue point was expected to be used to 
improve the photometric SED sampling of the targets. We 
explored the blue spectra in the 53-63 /im range (depend¬ 
ing on the source redshift) and found no clear signs of line 
emission in the data cubes, as expected. The significance 
of the continuum level (blue and red) was not high enough 
to permit a reliable photometric point for our SED fitting 
approach. 


4 RESULTS 

4.1 SED Properties 

Making use of the Herschel- PACS, Herschel- SPIRE, and 
the publicly available broadband measurements from WISE- 
22 /im (all targets detected) and IRAS (17 sources reliably 
detected at 60 /im), we fitte d the SEP of each galaxy follow¬ 
ing a similar approach as in llbar et al 1 (I 2013 I) . We fit a usual 
modified black body but forcing the SED shape to a power 
law in the high-frequency range of the spectra. The fits are 
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Figure 5. Dependency of the [Cll]/IR luminosity ratio as a function of dust temperature (top-left; colour-coded in mid-IR index), 
mid-IR index (top-right; colour-coded in /3) and dust emissivity /3 index (bottom; colour-coded in dust temperature). Data suggest that 
galaxies with colder dust temperatures tend to have higher [Cll]/IR ratios than hotter galaxies. 


aided by the known optical redshifts. We shift all SEDs to 
rest-frame frequencies, following (y — v 0 b s [1 +2]), and then 
fit the following parametrisation 

q ( \ _ a f MBB(z/) if v ^ v* 

b v (y) -AX| x {y/ Vt )- a mw-iR if v>Vt 

(1) 

where 


meaning, nevertheless it is useful to account for the the mid- 
IR part of the spectra which otherwise is underestimated by 
a simple modified black-body (see Appendix [A}. Examples 
of best fit SED are shown in Fig. [2] To measure the IR 
luminosity, L(8 — 1000/xm), we integrate the best fitted SED 
in rest-frame frequencies between — 0.3 THz (1000 /xm) 
and V 2 — 37.5 THz (8/xm), 


MBB(^) 


exp (fts) - 1 


( 2 ) 


The h and k parameters refer to Planck and Boltzmann 
constants, respectively. f3 is the usually called dust emis sivity 
index (e.g. ISeki &; Yamamotolll980l : iDunne et al.ll201lh - an 
averaged property of the dust grain emission over of whole 
galaxy. The parameter v * is obtained numerically at 


dlogio(MBB) 

diogioM 


(^*) — Q^mid-IR 


(3) 


which is simply used to match the slope of the modified black 
body function with the high-frequency power-law (roughly 
at ~ 100-200 jLim). This parameter does not have a physical 


L(8 — 1000 /im) = 47r D^(z) f S u dv (4) 

J 

Under this parametrisation, we obtain the dust tem¬ 
perature (Td), the dust emissivity index (/?), the mid-IR 
slope (<a m id-iR), and the normalisation which provides the 
total IR luminosity (rest-frame 8-1000/xm). Uncertainties 
for each parameter were obtained from a Monte-Carlo sim¬ 
ulation (lOOx), randomly varying the broadband photom¬ 
etry as appropriate for their measured uncertainties (Ta¬ 
ble [3]). The best-fit SEDs for all our target galaxies are 
shown in Fig. [2] We find that our sample have luminosi¬ 
ties of the order of 10 < logio(LiR,/L©) < 12, dust temper¬ 
atures of 20 < Ta/K < 55, mid-IR slopes in the range of 
1.5 < <a m id-iR < 3.2 (plus two outliers at ~4.5), and dust 
emissivity indices of 0.3 < < 2.5. 
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Figure 6. The figure shows the locus for our galaxies in a couple 
of colour-colour plots, exploiting Herschel and WISE photome¬ 
try. Given the optical extension of the targets, we consider ellip¬ 
tical aperture photometry for the WISE photometry. Broad-band 
luminosities have not been fc-corrected. Data points are colour- 
coded as a function of the [Cll]/IR luminosity ratio (see top-bar 
in left panel). We overplot also the expected c olour s for differ¬ 
ent SED templates taken from IPolletta et al.1 (l2007l Vs libraries 
(convolved with response filters) in five redshift ranges (z =0, 
0.05, 0.1, 0.15 and 0.2) indicated by small diamonds, where the 
biggest one is at z = 0. The figure indicates that galaxies with 
highest [Cii]/IR ratios are those with 0.5 < S'i2^m/5'22^m < 1.0, 
galaxies which could be explained by normal spiral ‘Hubble- class’ 
galaxies (e.g. SO, Sa, Sb, Sc), while those with lower values are 
better represented by star-burst templates such as Arp 220, M82, 
IRAS 22491, NGC6240, etc. 


4.2 The [C n] deficit 

Like many before us, we find that [C n] luminosity (including 
PDR and non-PDR components, which cannot be disentan¬ 
gled in this study) correlates strongly with Lir luminosity 
over approximately two orders of magnitude (Fig.0J). A sim¬ 
ple linear regression (excluding [C n] undetections) results in 
logio(L[ CII ]/L©) = (0.721±0.003) x log 10 (LiR/L©) + (0.53=L 
0.04) and the distribution has a Spearman’s rho rank corre¬ 
lation coefficient of 0.88 and a two-sided significance of its 
deviation from cero of 2 x 10 -7 . The slope of this correla¬ 
tion evidences the so-called [Cli] deficit, i.e. at higher IR 
luminosities the line to continuum ratio decreases. 


4.3 Correlations with dust properties 

The quality of our data permits us to explore dust temper¬ 
atures in great detail (uncertainties range from 1 to ~ 5 K 


based on Monte-Carlo simulations; see Table [3]). This is a 
major step forward compared to previous studies that only 
provided a proxy for the du st temperature, e.g. u sing the 
IRAS 100 /xm to 60 /xm ratio (|Malhotra et aljfeooil h Within 
our sample, most of the low-luminosity, low-redshift galax¬ 
ies with Lir < 10 11 L©, have isothermal dust temperatures 
of T d = 23 ± 1 K (Fig. 02) - typical of ‘normal’ H -ATLAS 
galaxies (e.g. ISmith et al.l 1201$ ). The galaxies with cooler 
dust (Td < 30K) tend to have ~ 2 x higher [Cll]/IR ra¬ 
tios than warmer galaxies. We find that ~ 90 per cent of 
the galaxies with Td < 30K have [Cn]/IR> 3 x 10 -3 . This 
fraction decreases to ~ 20 per cent for hotter galaxies. 

It is interesting to note the mild correlation where flat¬ 
ter mid-IR slopes tend to produce lower [C n]/IR ratios. This 
might be a manifestation for the presence of stronger radia¬ 
tion fields contributing to the 22 /xm flux density, or simply 
due to the poor mid-IR constraints on the SED-fitting ap¬ 
proach (see also g Qll . We find that the fitted dust emissivity 
index, /3, correlates with the [Cil]/IR ratio, where higher 
P values are found in galaxies with higher [Cll]/IR ratios 
(Fig. [5}. In Table [2] we explore various possible correlations 
for the [Cii]/IR ratio as a function of fitted and observed 
quantities. 


4.4 Correlations with WISE colours 

Fig. [ 6 ] shows a colour-colour diagram using WISE and Her¬ 
schel photometry for our sample. We find that low [Cii]/IR 
ratios (< 2.5 xlO -3 ) are almost ubiquitous among galaxies 
presenting a flux density ratio *Si 2 M m/< 0.5, i.e. a 
power-law with a^i < —1.14. Note that all of our galaxies 
got fitted <a m id-iR values which are steeper than —1.14 (see 
Fig. [5]), implying that the WISE colours do not reflect the 
same mid-IR slope. In fact, Fig. [ 6 ] indicates that the [C 11 ]/IR 
luminosity ratio depends strongly as a function of the 22 /xm 
flux density (see also Table [2]) . 

In Fig. [ 6 l we also show the expected WISE and Her- 
schel colours from different SED templates taken from 
IPolletta et al.l (|2007h and convolved with WISE transmis¬ 
sion filters in the redshift range of this study (z < 0 . 2 ). 
We find that mid-IR templates from star-bursting galax¬ 
ies (e.g. NGC6240) and moderately luminous AGNs (e.g. 
Sey 2 ) are well suited to explain low Si 2 ^m/S 22 ^m ratios (i.e. 
lower [Cn]/IR ratios). Higher [Cll]/IR ratios ratios tend 
to be better explained by spiral Hubble -type galaxies (S0- 
Sc). A close inspection to their mid-IR to far-IR SEDs, sug¬ 
gests three signatures controlling the S 12 /S 22^111 ratio; the 
strength of the PAH 7.7 /xm emission line system, the Sili¬ 
cate absorption band at 9.8 /xm and the slope of the spectra. 
Unfortunately, with the present data we are unable to distin¬ 
guish between these different features. We have performed 
a basic simulation taking into account the broadband fil¬ 
ters and artificial SEDs including different slopes, equivalent 
widths of the PAH emission and Si absorption band. We find 
that the mid-IR slope is probably the dominant component 
in varying the observed S 12 pm./S 22 pm. ratio. On the other 
hand, on average, galaxies with low (high) [Cn]/IR ratios 
have low (high ) 63 . 4 ^ 111 /S 4 . 6 ^m ~ 1.3 (1-45) ratios. Based on 
IPolletta et ahl ’s templates, we find that WISE colours from 
elliptical, or normal early-type galaxies, or powerful AGNs, 
lie outside Fig. [ 6 ] in agreement with the far-IR selection 
criteria which prefers star-forming dusty galaxies. 
























Table 2. The correlation of the L[cil]/-kiR luminosity ratio as a function of various different parameters. Linear regression values, where Lfcih/^IR = a-\-b X -0, are obtained by simple 
X 2 minimisation. We exclude the two [C il] undetections from this analysis. Last two columns show the Spearman’s rho (p) rank correlation and the two-sided significance of its deviation 
from cero (Q). We provide values for the best linear fit, only when a significant correlation is found Q < 5%. 


■0 

a 

b 

P 

Q{%) 

logioB'iRG'©)] 

-0.71214= 0.1021 

-0.17134= 0.0094 

-0.49 

1.1 

Tdust (K) 

-1.98034= 0.0180 

-0.01824= 0.0005 

-0.72 

0.0 

C*mid-IR 

-3.5550± 0.0256 

0.40824= 0.0105 

0.75 

0.0 

P 

-2.87054= 0.0141 

0.20844= 0.0094 

0.52 

0.7 

loglO [-^22/nm (erg/s/Hz)] 

4.93724= 0.2326 

-0.23884= 0.0074 

-0.76 

0.0 

logio R lOO/nm (erg/s/Hz)] 

1.4705=1= 0.3232 

-0.12314= 0.0098 

-0.42 

3.3 

logio R 160/u,m (erg/s/Hz)] 

- 

- 

-0.31 

12.8 

logio [^250/xm (erg/s/Hz)] 

- 

- 

-0.26 

19.2 

logio [-^350/nm (erg/s/Hz)] 

- 

- 

-0.26 

19.3 

logic R500/im (erg/s/Hz)] 

- 

- 

-0.25 

21.4 

^3.4/nm/^4.6/nm 

-3.4390=1= 0.0329 

0.62244= 0.0234 

0.62 

0.1 

l'12fim /-^22/u.m 

-2.94874= 0.0096 

0.84434= 0.0194 

0.87 

0.0 

-^22/itm/■I'lOO/um 

-2.29604= 0.0079 

-5.54474= 0.1360 

-0.82 

0.0 

r eff ,50% ( k P c ) 

- 

- 

0.19 

36.4 

r eff,90%( k P c ) 

- 

- 

-0.04 

82.8 

l°gl 0 [Si R (L o /kpc 2 )] 

-0.35704= 0.0792 

-0.24624= 0.0088 

-0.65 

0.0 


Table 3. Properties of the targets analysed in this work. Coll: H- ATLAS’s nickname; Col2 : GAMA’s redshift; Col3: Line peak flux; C 0 I 4 . : Intrinsic line FWHM width obtained after 
removing (in quadrature) the intrumental spectral resolution; Col5: Line flux density (\/27 t X P[cil] X fitted line width); C 0 I 6 : Line luminosity, where upper limits are 3 — a assuming 
an intrinsic 400kms _1 FWHM width; Col7: IR (8-1000pm) luminosity; C 0 I 8 : Dust temperature; Col9: Dust emissivity index; CollO: Mid-IR slope; Colli: Lfcih/^IR luminosity 
ratio, where upper limits use values from Column 6; Coll 2 : Effective radius at which 50% of the power is encircled using SDSS r-band imaging; CollS: Stellar mass; C 0 II 4 : Specific 
star-formation rate using the IR luminosity as proxy; Coll5: Metallicity measurements for star-forming galaxies; Coll 6 : BPT classification, 0=star-forming, l=composite and 2=AGN; 
Coll7: Optically based visual morphological classification, E=elliptical, S=spiral and I=irregular. 


Target 

z 

P [CII] 

FWHM [CII] 

S [CII] 

L [CII] 

I°gl0 [- L IR/- L Q] 

P dust 

P 

Q mid-IR 

[Cn]/IR 

r eff ,50% 
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sSFR IR 

12-|-logio (O/H) 

BPT 

Morph. 



(Jy) 

(km s — 1 ) 
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(xio -3 ) 

(kpc) 
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(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(7) 

(9) 

(10) 

(ii) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

G09.v2.102 

0.073 

1.904=0.06 

3564=22.88 

8494=46 

1.704=0.09 

11.034=0.01 

364= 2 

1.24=0.2 

2.104=0.06 

1.60 

3.5 

10.634=0.12 

0.30 

_ 

1 

ES 

G09.v2.103 

0.041 

4.454=0.06 

2204=0.10 

15094=31 

0.954=0.02 

10.284=0.03 

254= 2 

1.54=0.4 

3.004=0.18 

4.96 

5.0 

10.484=0.12 

0.08 

- 

1 

s 

G09.v2.107 

0.128 

4.504=0.09 

2274=21.31 

14674=62 

9.334=0.40 

11.704=0.01 

394= 1 

1.44=0.1 

2.504=0.05 

1.86 

4.2 

10.514=0.13 

1.84 

8.9234=0.011 

0 

I 

G09.v2.lll 

0.078 

3.394=0.07 

1924=18.32 

10524=42 

2.434=0.10 

11.154=0.01 

414= 2 

0.94=0.2 

2.304=0.06 

1.72 

5.1 

10.554=0.13 

0.48 

9.1474=0.013 

0 

ES 

G09.v2.117 

0.054 

2.994=0.06 

1684=12.24 

8984=31 

0.994=0.04 

10.424=0.04 

224= 2 

1.64=0.3 

2.504=0.18 

3.75 

7.6 

10.704=0.12 

0.06 

- 

- 

ES 

G09.v2.137 

0.044 

6.054=0.06 

2844=8.39 

23494=42 

1.714=0.03 

10.544=0.02 

344= 2 

1.34=0.2 

2.904=0.17 

4.90 

2.5 

10.134=0.13 

0.31 

9.0104=0.006 

0 

ES 

G09.v2.167 

0.078 

0.464=0.05 

- 

- 

< 1.25 

11.054=0.03 

504= 2 

0.74=0.1 

4.804=0.20 

<1.10 

2.5 

10.354=0.14 

0.61 

- 

2 

E 

G09.v2.170 

0.051 

4.524=0.07 

1954=12.97 

14394=38 

1.414=0.04 

10.334=0.03 

214= 2 

2.24=0.5 

2.504=0.13 

6.54 

4.4 

10.264=0.13 

0.14 

- 

- 

s 

G09.v2.175 

0.070 

3.034=0.06 

3124=17.73 

12354=46 

2.284=0.09 

11.184=0.01 

424= 2 

1.24=0.2 

1.904=0.05 

1.51 

3.4 

10.484=0.11 

0.60 

- 

1 

E 

G09.v2.232 

0.096 

2.154=0.07 

1844=0.01 

6484=32 

2.284=0.11 

10.934=0.04 

234= 3 

2.24=0.6 

2.704=0.20 

2.67 

7.5 

- 

- 

- 

1 

ES 

G09.v2.235 

0.027 

3.114=0.07 

1424=19.93 

9014=34 

0.244=0.01 

10.204=0.01 

504= 4 

0.34=0.2 

1.804=0.05 

1.53 

1.6 

9.944=0.11 

0.22 

8.6654=0.000 

0 

ES 

G09.v2.23 

0.033 

5.254=0.05 

3564=8.34 

23734=44 

0.974=0.02 

10.314=0.02 

244= 2 

1.64=0.4 

2.804=0.13 

4.72 

3.6 

10.554=0.11 

0.07 

8.9354=0.007 

0 

ES 

G09.v2.26 

0.182 

4.564=0.13 

1034=23.42 

10524=51 

13.814=0.68 

11.844=0.02 

244= 1 

2.14=0.2 

2.204=0.08 

2.02 

6.4 

11.254=0.12 

0.46 

9.1424=0.007 

0 

El 

G09.v2.299 

0.074 

2.444=0.06 

2584=25.34 

8804=48 

1.814=0.10 

10.984=0.01 

454= 3 

0.84=0.2 

2.404=0.08 

1.91 

2.3 

10.544=0.11 

0.33 

9.2144=0.012 

0 

E 

G09.v2.38 

0.059 

2.944=0.05 

5274=16.39 

17944=58 

2.304=0.08 

10.704=0.02 

224= 1 

2.04=0.3 

3.004=0.17 

4.55 

4.6 

10.894=0.11 

0.08 

9.1154=0.648 

0 

ES 

G09.v2.42 

0.055 

3.894=0.07 

2154=17.74 

12914=43 

1.454=0.05 

11.094=0.01 

384= 1 

1.34=0.1 

2.304=0.04 

1.17 

2.8 

10.464=0.12 

0.51 

- 

1 

E 

G09.v2.45 

0.051 

6.104=0.06 

2664=7.28 

22744=41 

2.184=0.04 

10.714=0.01 

314= 1 

1.34=0.2 

3.104=0.13 

4.27 

4.0 

10.284=0.13 

0.32 

8.9014=0.009 

0 

s 

G09.v2.48 

0.072 

7.174=0.07 

2234=7.67 

23984=43 

4.694=0.09 

11.274=0.01 

364= 1 

1.44=0.1 

2.104=0.03 

2.51 

2.4 

10.584=0.13 

0.59 

- 

1 

El 

G09.v2.52 

0.026 

4.354=0.06 

1684=11.33 

13314=33 

0.334=0.01 

10.254=0.01 

274= 1 

1.84=0.3 

1.904=0.05 

1.81 

2.7 

10.294=0.11 

0.11 

- 

2 

ES 

G09.v2.55 

0.054 

4.284=0.07 

2964=15.46 

16974=53 

1.864=0.06 

10.574=0.02 

234= 1 

2.04=0.3 

2.904=0.11 

5.05 

3.3 

10.814=0.11 

0.07 

- 

1 

ES 

G09.v2.58 

0.052 

5.414=0.06 

2234=9.54 

18324=39 

1.864=0.04 

10.694=0.02 

254= 2 

2.04=0.4 

2.604=0.10 

3.81 

3.4 

10.604=0.12 

0.15 

9.0864=0.047 

0 

ES 

G09.v2.60 

0.060 

3.524=0.05 

4824=13.53 

19974=55 

2.664=0.07 

10.724=0.02 

234= 2 

2.04=0.4 

2.904=0.17 

5.03 

5.0 

10.634=0.12 

0.15 

- 

1 

s 

G09.v2.66 

0.031 

6.734=0.07 

1774=8.13 

20904=35 

0.734=0.01 

10.134=0.03 
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4.5 Correlations with optical properties 


0.010 


In this section, we explore the optical properties of the tar¬ 
gets. First, we look for morphological features. We down¬ 
loaded the corresponding fits files r-band images from the 
SDSS archival (see Fig. [3]). We masked all sources which 
are not contributing to the far-IR emission, and created the 
encircled energy fraction as a function of a radius centred at 
the peak source position. We measured the effective radius 
at which 50 % (and 90 %) of the power is encircled, and then 
transformed this projected size into a physical scale at the 
given redshift of the galaxy. These scales are shown in Fig.0 
as a function of the [Cii]/IR ratio. No clear correlation is 
found between these two observables (Table 0. 

We then performed a basic visual inspection using the 
r-band SDSS imaging to classify these galaxies in three pop¬ 
ulations: ellipticals (‘E’), spirals (‘S’) and irregulars (‘I’). 
Some galaxies show morphologies which are not possible to 
classify in a single population. In these cases we use a com¬ 
bination of letters, e.g. ‘ES’ or ‘EF, to refer to a prominent 
bulge with a disk or an irregular morphology, respectively. 
Under this classification, in Fig. [7] we show that the [C n]/IR 
ratio is preferentially higher in galaxies presenting a promi¬ 
nent disk compared to those which do not present disky 
morphologies. 

In order to provide an estimate for the IR surface 
brightness of the galaxies, we consider the optical r-band 
radius as a proxy for the actual IR radius, so we define 
Eir = LiR/( 27 rrg ff 50% ). We reckon this is a strong assump¬ 
tion so this value should be interpreted with caution. As 
expected, we find that the previous morphological classifi¬ 
cation is correlated to the surface brightness. Indeed, those 
with low Eir values tend to be those presenting disk-like 
morphologies, while those which are classified as pure ellip¬ 
tical tend to be those with highest brightness in the sample. 

With the aim to ident i fy po ssible [Cn] self-absorption 
(e.g. following iGerin et al.l 12015 1. we have looked for the 
properties of edge-on galaxies within the spiral population 
(see Fig. [3]). The four identified edge-on galaxies tend to 
have a well defined high [Cli]/IR ratios of 3.5 x 10 -3 , and 
do not deviate from the rest of the ‘spiral-only’ population. 

To characterise the nature of our galaxies, we have 
utilised the SDSS and GAMA spectra to locate them 
in they “B aldwin, Phillips & Terlevich” (BPT) diagram 
(Bal dwin et al.lfl98lh . Emission line strengths for H<a, H/3, 
[N n] A6583 and [O ill] A5007 are extract ed from spectra usin g 
the GANDALF pipeline (see details in lHopkins et al.ll2013h . 
Only two galaxies do not have all four lines detected. We 
show these data in Fig. [8j where our sample is over-plotted 
on top of the whole spectroscopic SDSS population. We find 
that our sample locates all over the star-forming, AGN and 
composite (a mix of both) regions. The main selection cri¬ 
terion, *Si 60 Mm > 150 mJy, therefore introduces a significant 
number of low-luminosity AGNs in the sample. This per¬ 
mits to test if AGN activity plays an important role in the 
[Cll]/IR ratio. Note that none of our galaxies falls close to 
the peak of the distribution of star-forming SDSS galaxies. 
We show that under the BPT classification, our sample of 
star-forming and composite galaxies do not show significant 
differences in terms of [Cll]/IR ratios (Fig. [8]). We note, 
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Figure 7. Top: [C n] /IR luminosity ratio as a function of effective 
radius in SDSS r-band (encircled fraction at 50% in black and 
90% in red). Bottom: The variation of the [Cll]/IR luminosity 
ratio as a function of surface brightness (Sir = Fir/(27t r^ ff 50% ); 
here r e ff 50 % is based on the optical r-band size and Eir has not 
been corrected by inclination), colour coded as a function of a 
simply by-eye optical morphological classification; spirals (‘S’), 
ellipticals (‘E’), irregulars (‘I’; including mergers), or composite 
morphologies. This figure clearly shows that higher [C 11 ]/IR ratios 
are preferentially seen in galaxies presenting extended disks. 


however, that the uncertainties in the line flux ratios might 
be blurring any possible correlation. 

The stella r masses (MG for a ll galaxies are calculated 
as described in lTavlor et al.l (l201lh . using GAMA catalogue 
version v08. The s tellar mass estimates we re derived by 
fitting their SEDs ([Bruzual Sz Charlotl l2003h to the SDSS 
ugriz imaging - data which ha ve been re-processed by the 
GAMA team ([Hill et al. 2011). Th e dus t obscur ation law 
applie d was that of ICalzetti et al.l ()2QQQh , and a IChabrieJ 
12003 1 IMF was assumed. The stellar masses are determined 
by integrals that are weighted to the probability of each 
SED fit. This has been performed to all galaxies regard¬ 
less of the nature (star-forming/AGN) defined by the BPT 
diagram. We have not applied any conversion factor to con¬ 
vert from Chabrier to Kroupa IMF because the variation 
is negligible compared to the measured errors. iTavlor et al.l 
demonstrate that the relation between (g-i) and M*/L of¬ 
fers a simple indicator of the stellar masses. To account for 
aperture effects, a correction based on the Sersic fit to the 
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Figure 8. Top: The nature of our targets as shown by the 
BPT diagram (see § 14.51) . Nine galaxies present line ratios with 
large error bars (0.1 dex). Ga laxies to the left of the solid line 
(taken from iKauffmann et ahl 120031 ) are considered star-forming 
( blue crosses); thos e to the right of the dashed line (taken from 
iKewlev et al .11200lh are considered AGNs (red crosses); and those 
between the solid and dashed lines are defined as composites 
(green crosses). The background contours and grey scale shows 
the distribution for the whole spectroscopic SDSS sample. Bot¬ 
tom: The [C ll]/IR luminosity ratio as a function of IR luminosity, 
colour-coded by the BPT classification (same colours as above). 
The present data show that in this diagram the different popula¬ 
tions are indistinguishable from each other. 

surf ace brightness profiles is applied to t he stellar masses 
(see iTavlor et al.l 1201 ll : I Kelvin et al.l l2012h . 

We use th e SF R derived fro m the bolometric IR emis¬ 
sion, following iKennicii tt] ill 99811 . In combination with the 
stellar mass estimate we obtain the specific-SFRs (sSFR = 
SFRir/M*). We find the sSFR seems to anti-correlate with 
the [Cii]/IR luminosity ratios (see Fig. [9}, hence those 
galaxies presenting lower [Cll]/IR ratios are preferentially 
passing through more violent bursts of star-format ion. In 
order to remove a possi b le dep endency on redshift and foil¬ 
wing iDiaz^Santos^l^l] (|2013h ’s work, we have divided our 
sSFR esti mates by the one e xpected at the ‘main-sequence’ 
defined bv lElbaz et al.l (|2Q1 lh (at the given redshift). We find 
that our sources are centred at the main sequence value but 
cover a wide range of two orders of mag nitude around it (s ee 
Fig. [9]). If we compare our sample to iDfaz-Santos et alJ ’s. 


our sample is more representative of ‘normal’ star-forming 
galaxies than theirs. 

For those galaxies of which are defined as star-forming 
in the BPT diagram (Fig. [8}, we measure relia b le me tallic- 
ities using the o3n2 index of IPettini &; Pag e] (I2QQ4I ) , and 
conve rting those v alues to the calibr a tion o f Tremonti et al .1 
(l20Q4h , following lLara-Lopez et al . (|2Q13h . We find that 
these star-forming galaxies present metallicities between 
12 + logio(0/H) = 8.7 — 9.2. For comparison, Solar 

metallicity is 8.91. In low metallicity environments (12 + 
logio(0/H) < 8.1), the IR emission is expected to drop but 
the [C n] emission would remain almost invariant, therefore 
the [Cn]/IR ratio is exp ected to be higher in low metallic¬ 
ity en vironments (see also lRubin et a l. 2009; D e Looze et al.l 
l2014h . Nevertheless, our sample does not probe these lower 
metallicities (mainly due to dust selection criterion) and so 
we cannot investigate correlation of [C n] with metallicities. 


5 DISCUSSION 

For many years it has been a great challenge to explain why 
there is a decrease in the [Cn]/IR luminosity ratio towards 
bright (Lir > lO 11 !/©) IR luminosities. It is well established 
that this ratio has an intimate dependency on the strength 
of the radiation field (<Go>) and the density (n e ) of the 
ISM. As explained in §[l] the [C n] emission comes from a 
whole range of different ISM states; ionised, atomic and/or 
molecular. Unfortunately, given the low resolution of our 
Herschel observations, these different components are im¬ 
possible to disentangle directly from imaging. Fibre-based 
optical spectroscopy cannot spatially separate these compo¬ 
nents neither. 

First of all, we look for possible biases introduced by 
the selection criteria used to construct our galaxy sample. 
We identify that one of the most important criterion is the 
flux density Si6o^ m >150 mJy threshold, which introduces 
a strong selection effect on the redshift-luminosity distribu¬ 
tion of the targets (see Fig. [lj . Our targets span a redshift 
range between z — 0.02 and 0.2, so approximately 2.2 Gyr 
of cosmic time. Based on the measured sSFR, the galax¬ 
ies in our sample double their stellar masses in scales of 
sSFR -1 « 0.3 — 20Gyr (median 4Gyr; see Fig. [9}, values 
which are in most cases longer than 2.2 Gyr. In this work, 
we consider that for most of our targets the evolution is 
short in comparison to the redshift slide, i.e. a galaxy at 
fixed luminosity will be behaving the same at z = 0.02 as 
at z — 0.2. Only a fourth of the targets double their stel¬ 
lar masses in scales which are shorter than 2.2 Gyr (those 
ongoing more violent star-bursts). This assumption is nec¬ 
essary to alleviate the luminosity-redshift dependency seen 
in Fig. [T] In te rms of the ISM evolutio n using large sam¬ 
ples of galaxies, lLara-Lopez et al ] ll200flh finds no variation 
of the metallicity properties in this redshift range, neverthe¬ 
less the dust mass density evolves strongly as a function of 
reds hift, incrementing by a factor of two from z = 0.02 to 
0.2 (| Dunne et al.ll201ff ). For the purposes of this work, we 
assume the flux density threshold at S mo^m > 150 mJy does 
not introduce biases on the results. 

Significant differences are found between galaxies pre¬ 
senting high and low [Cli]/IR luminosity ratios (see Ta- 
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Figure 9. Top: The [Cll]/IR luminosity ratio as a function of 
the sSFR, derived using the observed IR luminosity and stel¬ 
lar masses from SEDs fits to SDSS ugriz photometry (see § 14.51 
for details). Colour-coded we show metallicity measurements (us¬ 
ing [TremontU^d] EoCmI’s relation) for only those sources which 
are identified as star-forming in the BPT diagram (Fig. El). Bot¬ 
tom: The same figure but dividing the sSF R values by the red - 
shift dependent ‘main-sequence’ defined by lElbaz et al ] <201 ill . 
In light grey aster i sks w e show the galaxy sample presented by 
IDfaz-Santos et a hl (bo d). The dashed line shows the correlation 
defined bv iDfa^Santos et al.l (see their Eqn. 5 but converting 
L m = 1-75 x Lfir), including a =L0.26dex range of dispersion 
(dotted lines). Note how our sample is more representative of 
‘normal’ galaxies (wi th respect to the ‘main-se quence’) compared 
to that presented bv I Diaz-Sant os et al ] d2013l ). 


ble [2]) - understood in this study as greater or lower than 
~2.5 x 1CU 3 . We find that galaxies with high ratios: 

• have cold dust temperatures, preferentially lower than 
30 K (see Fig. El). This evidence indicates that these galaxies 
present a prominent IR component coming from extended 
ISM regions rather than compact ones located at the vicinity 
of powerful star-forming nuclear regions; 

• have a high WISE flux density ratio of 0.5 < 

iSi 2 /xm /S 22 iim <1.0 (see Fig.EJ). This range of mid-IR ratios 
tend to be associated to normal Hubble -type spiral galaxies, 
a ratio induced by a combination of spectral features within 
the broad WISE 12 /im and 22 /im filters: prominent PAH 
emission lines (not suppressed by strong radiation fields), 
a weak 10 /im Silicate absorption band (indicating moder¬ 


ate extinction levels), and a ~ 22 /im spectra which does not 
seem to be dominated by powerful hot dust emission, e.g. a 
young and violent star-burst or an AGN torus; 

• have preferentially lower surface brightness (Eir < 
10 9 L©/kpc 2 ) as shown in Fig. [7| (see significance in Table[2|). 
We stress this estimate uses the effective r-band radius as 
a proxy for the IR extension. If the mean free path of the 
far-UV photons is large, then the strength of the radiation 
fields would be directly proportional to Eir (|Wolfire et al.l 
Il990l ). If this is the case, this suggests that galaxies with 
weaker radiation fields produce higher [Cll]/IR luminosity 
ratios. Nevertheless, given the lack of correlation between 
[Cli]/IR and r e ff (see Table [2]), this behaviour could be a 
manifestation of the correlation with Lir instead; 

• have preferentially disk-like morphologies. Those galax¬ 
ies which are classified as spirals (‘S’), without prominent 
bulges, are those with the highest [Cll]/IR ratios 4 x 
10~ 3 ; see Fig. 0 , similar to that found bv lDfaz-Santos et al.l 
(120131 ) while looking at the pure star-forming LIRGs in their 
sample. These evidences suggest that galaxies with high 
[Cii]/IR ratios evolve quiescently rather than triggered by 
a major merger event with a subsequent powerful nuclear 
star-burst; 

• present a wide range of sSFR^ 5x 10“ 11 — 3x 1CU 9 yr -1 
(IR-based), hence the mechanism controlling the [C 11 ] emis¬ 
sion does not seem to relate to the efficiency of convert¬ 
ing gas into stars. At constant SFR, these galaxies range 
from 3.3 to 20 Gyr to double their stellar masses. These re¬ 
sults do not agree with the [C 11 I /IR versus sSFR/sSFRm s 
correlation (see Fig. [9]) found bv IDfaz-Santos et al.l (120131 ). 
probably because we observe more ‘nor mal’ galaxies (rel - 
ative to the ‘main-sequence’ defined by lElbaz et al.l |201l] ) 
than their local LIRGs sample, and also given by the fact 
that our sample does not include a significant number of 
galaxies with low [Cll]/IR ratios that permit to evaluate 
the correlation at higher sSFR levels. Our data m ight sug¬ 
gest an evidence for a plateau in lPfaz-Santos et ahl ’s correla¬ 
tion at sSFR/sSFRms < 1, probably induced by a different 
and more inefficien t star-format ion mechanism controlling 
the [ C 11 I/IR ratio (IPaddi et al.l l2010l : iGracia-Carpio et al.l 

ImT 


5.1 The strength of the radiation field 

Modelling the [C 11 ] emission as coming from PPRs, we sug¬ 
gest that one of the main parameters responsible in control¬ 
ling the [C ll]/IR ratio is the strength of the far-UV radiation 
field (<Go>). This is supported by the significant correla¬ 
tion found between dust temperatures, and Eir, with the 
L[ch]/Tir luminosity ratio. Higher dust temperatures sug¬ 
gest higher radiation fields generated by higher SFR surface 
densities, which might create large ionised complexes, espe¬ 
cially expected in those galaxies with high IR luminosities 
11 < logio(TiR/L©) < 12. These more extreme conditions 
could easily change the dominant ISM state responsible fo r 
the bulk of the [Cn] emission (e.g. IPfaz-Santos et al.ll2013l ). 

Actually, the far-UV radiation field (produced by O 
and B stars) is one of the main contributors to the heat¬ 
ing of the gas via the photo-electric effect on dust grains. In 
the case of soft radiation fields, t he ejection rate of photo¬ 
electrons from dust decreases fe.g. lSoaans et "ahll 19941 ), while 
in the case of strong radiation fields, the dust grains be- 

















































































come positively charged, increasing the potential well that 
the photo-electrons need to overcome, and thus reduc¬ 
ing the i nput energy transferred to the gas by photo¬ 
electrons (iTielens Sz Hollenbachl[l985l : iMalhotra et al.)[l997l . 
l200ll : [Luhman et al.l 2003h . 

We suggest that the [Cii]/IR ratio is controlled by the 
strength of the far-UV radiation fields, hence the decrement 
of the [C 11 ] line with respect to IR emission is most probably 
due to an increment of positively charged dust grains (higher 
dust temperatures), which reduces the efficiency of the far- 
UV radiation field in tra nsferring energy into the gas. 

iNegishi et. al] ll200lh found that the <Go>/n e (where 
n e is the density of electrons of the ISM) ratio does not drive 
the [Cii]/IR ratio but they suggest that high gas densities 
play an important role in controlling the [C 11 ] emission. Mor¬ 
phologically speaking, our analysis shows that galaxies with 
low [C 11 ] /IR ratios tend to have prominent bulges in nuclear 
regions, i.e. probably suggesting that gas density plays an 
important role in the [Cll]/IR ratio. However, we were un¬ 
able to identify a clear correlation between [Cii]/IR ratio 
and sSFR, indicating that more efficient SFR in compact 
regions is probably not controlling the [Cll]/IR ratio (at 
least at the parameter space explored by this work). Un¬ 
fortunately, with the available data presented in this work, 
we are unable to separate the intimate relation between the 
strength of the radiation field and the density of the ISM. To 
separate both parameters we require [C 11 ] together with an¬ 
other emission line, such as the fine transitions of [N 11 ], [O 1 ] 
and [C 1 ] or rotational transitions of CO to p roperly deter¬ 
mine the physical conditions of the ISM (e.g. IWoffire et al.l 
Il989l : iHailev-Dunsheath et al.l 2010 ). 

5.2 The ISM origin of the [C 11 ] emission 

We show that galaxies presenting high [Cii]/IR ratios have 
relatively cold dust temperatures, have dominant disk-like 
morphologies, and low surface brightness, evidences that in¬ 
dicate relativel y weaker far-UV rad iation fields. As previ¬ 
ously shown bv IPineda et al.l (2 014b . the origin of the [Cll] 
emission in the Milky Way (a ‘normal’ spiral galaxy) is not 
only from cold PDRs, but includes also contributions of the 
same order from ionised gas, diffuse atomic gas and CO-dark 
H 2 . It is expected that weaker radiation fields would imply 
smaller complexes of ionised gas. We argue that sources with 
high [C ll]/IR ratios might not only emit their [C 11 ] luminos¬ 
ity from cold PDRs but also from the diffuse and extended 
atomic ISM phase. Without a H 2 tracer (e.g. CO lines), we 
are unable to prove this statement, although it points out 
to the difficulties in understanding the origin of the [C 11 ] 
emission with a single far-IR line detection. 

5.3 Old stellar populations contributing to the 
cold IR emission 

The [C 11 ] and IR luminosities are intimately related to the 
star-formation process. In this section we explore if the 
‘[C 11 ]-deficit’ could be a manifestation of an inclusion of an 
IR emitting component which is not related to the star- 
formation, but to old stellar populations. Actually, the IR 
SED component coming from dust heated by old stellar pop¬ 
ulations is cold, diffuse and is predominantly emitted at long 


> 200 fim wavelengths. We have explored the possibility that 
the 500 jam luminosity could correlate with the [C ll]/IR ra¬ 
tio, although no clear trends are observed (see Table [2]). If 
a prominent cirrus-like emission is present in these galax¬ 
ies, the Rayleigh-Jeans would tend to have flatter spectra 
(between 250 and 500 /mi), hence would tend to bias the 
fitted dust emissivity index (see also Appendix [Aj - this pa¬ 
rameter basically controls the slope in the Rayleigh-Jeans 
regime. With the available data it is not possible to distin¬ 
guish different physical dust properties from a strong cirrus 
component. Nevertheless, it is worth mentioning that the 
anti-correlation found between best fit ft and the [Cii]/IR 
luminosity ratio (Fig. [5]) suggests that it is unlikely that the 
bolometric IR luminosity is dominated by a cold cirrus com¬ 
ponent, hence responsible for the decrement seen in [C 11 ] /IR 
as a function of IR luminosity. This is in agreement with the 
selection criterion SiQo^m > 150 mJy which prefers galaxies 
dominated by star-forming heating. 


5.4 Self absorption or optically thick 

Extraordinarily large column densities are required to make 
the [C 11 ] emission optically thick. Self absorption has 
been employed to explain the [C III emission from AGN- 
dominated systems, like Mrk231 (iFischer et al.l 20101 where 
Hi column densities could be higher than 10 22 cm -2 . For 
most star-burst galaxies this effect is expected to be small, 
especially at the luminosity range explored by this wor k 
10_< logio(T ir/-£/©) < 12 . On the other hand. [Germ et all 
(l2015h have recently shown that on the plane of the Milky 
Way the presence of foreground absorption may completely 
cancel the emission from a background far-IR emitter in 
medium spectral resolution data, suggesting that spectra 
should be taken at high spectral resolution, e.g. using HIFI 
(Ide Graauw et aflEncn rather than PACS, to interpret cor¬ 
rectly the [Cn] emission, therefore the [Cil]/IR ratio. 

In order to explore this idea, we have identified all four 
edge-on galaxies in our sample as possible optically thick 
[Cll] candidates (along the line of sight). We find that all 
four galaxies are at the high end of the [C 11 ] /IR luminosity 
ratio distribution. This result suggests that the [Cll] emis¬ 
sion might not be self-absorbed, at least by the disk, where 
most of the Hll regions and surrounding PDRs are placed. 

Supposing that the whole [Cll] emission of the disk 
(mostly PDR-related) is absorbed, then the observed high 
[Cii]/IR ratios should come from a diffuse [Cll] component 
preferentially located above/below the disk. Note also that 
these edge-on galaxies tend to show line FWHMs which are 
> 200 kms -1 (see Table [3l) . helping the line emission to es¬ 
cape from the disk. 

5.5 IR emission contaminated by an AGN 

It has been previously shown that sources harbouring 
an AGN have lower [C 11 ] to IR luminosity ratios (e.g. 
ISargsvan et al.l l2012h . In Fig. [8] we have classified star¬ 
forming galaxies from AGNs using the BPT diagram. Error 
bars are large, although we find that star-forming, compos¬ 
ite and AGN populations are indistinguishable in terms of 
[C 11 ] /IR luminosity ratios, possibly suggesting that the pres¬ 
ence of an AGN might be playing a local but not a global 
































role on the [Cll]/IR luminosity ratios - at least at the low 
A GN luminosities pre s ented in this work. This is supported 
by IDfaz-Santos et al ] (l2013h as they found that AGNs se¬ 
lected by a simple PAH equivalent width threshold do not 
modify significantly the [Cii]/IR luminosity ratio. 


6 CONCLUSION 

We have used recent PACS spectroscopic [C il] observations 
to describe its relation to the IR luminosity in a sample of 
28 galaxies selected from the if-ATLAS survey. This sam¬ 
ple has high-quality IR photometry from WISE, IRAS and 
Herschel , with the addition of unambiguous photometry and 
spectroscopy from the SDSS and GAMA surveys. 

In summary, after an exploration over a wide multi¬ 
wavelength parameter space, we have identified the following 
correlations. We find that galaxies with high L[ C ii]/Tir > 
2.5 x 10 -3 luminosity ratios tend to: have Lir < 10 11 L 0 , 
dust temperatures lower than 30 K, high WISE colours in 
the range 0.5 < S'i2 At m/*S22^m < 1-0, present disk-like mor¬ 
phologies, have low surface brightness Sir ~ 10 8-9 L 0 /kpc 2 
(using the r-band effective radius), and got a range of star- 
formation rate efficiencies (sSFR~0.05 — 3Gyr _1 ). 

Assuming that the physical properties of star-forming 
galaxies, at fixed luminosity, are the same at z = 0.02 and 
0.2 (the range of redshift of our galaxy sample), and based 
on the correlations found between the [Cii]/IR luminosity 
ratio and the dust temperature (and Sir), we conclude that 
the most probable parameter controlling the [Cii]/IR lumi¬ 
nosity ratio is the strength of the radiation field (averaged 
over the entire galaxy) - probably inducing an increment 
of the positive charge of dust grains that has an effect in 
the effective energy deposited by the radiation held in the 
electrons extracted from dust grains. 

The lack of correlation between galaxies with high 
[Cii]/IR luminosity ratios (> 3 x 10 -3 ) and sSFR values 
suggests that the efficiency to convert gas into stars (e.g. in 
high density environments) is not playing a dominant role in 
the line t o continuum behaviour, contrary to the correlation 
found bv lDfaz-Santos et al.l |2Q13 b We are probably observ¬ 
ing a plateau for the correlation at lower SFR/SFRms < 1 
ratios, maybe product of a different star-forming mechanism 
that controls the [Cll]/IR ratio. 

We find that the [C il] deficit is unlikely to be a man¬ 
ifestation of optically thick [C n] emission, as evidenced by 
the high [Cii]/IR luminosity ratios found in edge-on spiral 
galaxies. On the other hand, the analysis we performed to 
characterise the nature of our galaxies, using the BPT dia¬ 
gram, is not conclusive as all star-forming galaxies, AGNs 
and composite populations do not clearly distinguish from 
each other in terms of the [C il] /IR luminosity ratio. We con¬ 
clude that at least at the AGN luminosities shown by our 
sample, AGN activity does not seem to play a dominant role 
in the [C n] deficit. 
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APPENDIX A: 

TESTING FITTED PARAMETERS 

As described in §[4j our SED fitting approach includes a 
non-standard method. The inclusion of the power-law in the 
mid-IR forces the slope of the modified black-body emission 
at ~ 100-200 /im, possibly introducing a bias on the derived 
dust temperature or dust emissivity index. For this reason 
we repeated the SED fitting approach excluding the mid- 
IR slope, leaving just the modified black body (MBB from 
Eqn[2]) component. For these purposes we just use the Her- 
schel photometry, i.e. the 100, 160, 250, 350 and 500 /im data 
points. We also restricted, between 1.5 and 2.5, the range of 
possible values for the dust emissivity index. 

The first thing to note is that these new fits are unable 
to describe the high-frequency part of the spectra, hence 
we cannot use them to get bolometric IR measurements. 
These new fits, however, show a clear difference between 
derived parameters. On average, we find that with this new 
SED-htting method, the Td us t decreased by ~6K while ft 
increased by ~ 0.5 - parameters which are well know n to be 
correlated ie.g. lShettv et al.ll2009l : ISmith et al.ll2013h . These 
results are shown in Fig. iah This clearly demonstrate that 
converting from fitted parameters to ‘physical’ parameters 
should be taken with great caution. 

Even though a significant difference is seen between de¬ 
rived parameters coming from these two different SED fit¬ 
ting approaches, the previous trend seen in [Cii]/IR lumi- 
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Figure Al. Top: A comparison between derived parameters, 
^dust and P, using two different SED fitting approaches. In grey, 
we show the values obtained following the method described in 
§IH In red, we present the values obtained by fitting a modified 
black body emission using only the Herschel 100-500 p,m pho¬ 
tometry (as usually performed in previous Herschel-b&sed stud¬ 
ies). Bottom: The [C ll]/IR luminosity ratio as a function of fitted 
dust temperature. Colours are the same as in the top figure. This 
comparison shows that the dependency for [C il] /IR as a function 
of Td us t (seen in Fig.[5j remains with this different SED-fitting 
approach. 


nosity ratio as a function of dust temperature remains (see 
Fig. IA1D . This confirms, again, that high [Cii]/IR ratios are 
associated to galaxies dominated by cold dust emission. 





























